Frustrated antiferromagnets at high fields: the Bose-Einstein condensation in 

degenerate spectra 
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Quantum phase transition at the saturation field is studied for a class of frustrated quantum anti- 
ferromagnets. The considered models include (i) the J1-J2 frustrated square-lattice antiferromagnet 
with J2 = |Ji and (ii) the nearest-neighbor Heisenberg antiferromagnet on a face centered cubic 
lattice. In the fully saturated phase the magnon spectra for the two models have lines of degenerate 
minima. Transition into partially magnetized state is treated via a mapping to a dilute gas of hard 
core bosons and by complementary spin-wave calculations. Momentum dependence of the exact four- 
point boson vertex removes the degeneracy of the single-particle excitation spectra and selects the 
ordering wave- vectors at (n, tt) and (tt, 0, 0) for the two models. The asymptotic behavior of the mag- 
netization curve differs significantly from that of conventional antiferromagnet in d-spatial dimen- 



sions. We predict a unique form for the magnetization curve AM 
where is a distance from the quantum critical point. 
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Heisenberg antiferromagnets exhibit a quantum phase 
transition between a fully polarized state and a state with 
partial magnetization at the saturation field H c . A par- 
tially magnetized state typically breaks spin-rotational 
symmetry about the field direction and has a long-range 
order of transverse spin components. Such a zero tem- 
perature transition belongs to the XY universality class. 
Various properties of quantum antiferromagnets in the 
vicinity of this quantum critical point are well understood 
in terms of the Bose condensation of magnons below H c 
0, Q| ■ The antiferromagnetic wave- vector corresponds to 
the minimum in the magnon spectrum, whereas an arbi- 
trary phase of the condensate describes sublattice orien- 
tation in the plane perpendicular to the field. 

The above simple picture fails, however, for so-called 
frustrated antiferromagnets, which have degenerate clas- 
sical ground states between zero and saturation fields 
0, . The excitation spectra of frustrated antiferromag- 
nets above H c are quite unusual. The magnon dispersion 
e q has a continuous set of energy minima given by a d'- 
dimensional hypersurface embedded in the d-dimensional 
reciprocal space. The problem of the Bose condensation 
of particles with such degenerate spectra remains unex- 
ploited to a large extent. Numerical investigations 0,0 
of a specific model (model A below) show a strong sin- 
gularity of the magnetization near the saturation field, 
which has not been interpreted so far. In the present 
work we investigate two problems: (i) how the quantum 
fluctuations remove degeneracy of the excitation spectra 
and select a certain ordering wave- vector and (ii) how the 
asymptotic behavior of the magnetization curve close to 
the quantum critical point is modified by a presence of 
large phase space of soft-mode fluctuations. Specifically, 
we consider two models with d! = 1 in d = 2 and d = 3 
spatial dimensions: (model A) the 3\-3i antiferromag- 



netic Heisenberg model on a square lattice at its criti- 
cal point J2 = \ J\, which has lines of minima at (jr, q) 
and (q, tt) and (model B) the nearest-neighbor Heisen- 
berg model on a face centered cubic lattice (fee) with 
minima at (tt, q, 0) and equivalent lines. 

Convenient approach to deal with quantum antiferro- 
magnets near the saturation field is to employ a hard- 
core boson representation of spin- 1/2 operators: Sf — 
h - b\b t , S+ = b,, and Sr = b\ with b\b t = or 1 
0. The hard -core constraint is imposed by an on-site 
repulsion U -^00. The Heisenberg spin Hamiltonian 
H = J2u j) JijSiSj ~ I s ! then, transformed to 
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where e q = i(7 q -7„ 
transform of the exchange interaction Jy , /i = H c — H is 
a boson chemical potential, and H c = ^(70 — 7min) is the 
saturation field. A bare four-point boson vertex is given 
by = C^ + 7q- The exact ground state of the Hamilto- 
nian is the boson vacuum for ^ < (H > H c ), it cor- 
responds to a ferromagnetic alignment of spins. In d > 2 
the ground state with finite boson density (b}^) ^ at 
fi > should be, generally, a superfluid state: (6 q ) ^ 
for a particular wave- vector q = Q, such that £q — fx 
vanishes at the saturation field. 

In the above two models the bare spectrum e q has 
lines of degenerate minima (see below) and the conden- 
sate wave- vector remains undetermined. The degeneracy 
cannot be lifted by the bare interaction V q as it depends 
on momentum transfer only. This is a manifestation of 
an infinite classical degeneracy of the ground state of the 
original spin problem below the saturation field. Such 
an 'accidental' degeneracy can be removed by quantum 



7 q = V ■ Jije l<Vi i is the Fourier 
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fluctuations jg ■ The number of bosons vanishes at the 
quantum critical point and hence the system becomes 
dilute in the limit /i — » + . In this case the leading or- 
der corrections to the scattering vertex are generated by 
the processes in the particle-particle channel. According 
to the standard procedure the bare interaction is, then, 
replaced by a solution of the integral equation for the 
scattering amplitude T |<| . We show that T does depend 
on the incoming momenta and thus lifts the degeneracy 
of the original problem. The ordering wave- vector corre- 
sponds to the smallest T, i.e., condensate occurs at the 
momentum at which bosons less interact with each other 
in order to minimize their repulsion. 

The Bethe-Salpeter equation for the scattering func- 
tion with zero total frequency reads as 



r q (k,k') = v q -i^v q _ p 



r P (k,k') 

e k+p + Ek'-r 



(2) 



and is graphically represented in Fig. ^ In the limit 
U — > oo, Eq. (01 is reduced to a system 



r q (k,k') = 7q + (r)-^7q-p 



r P (k,k) 
v^ q - p ek+P + e k'-p 



iy r p (k,k) _ 1 

N ^ £ k +p + £k'-p 



(3) 



where (r) = (1/N) E q r q( k > k and the identity (7) = 
Ju = has been used. By expanding r q (k, k') in lattice 
harmonics of the wave-vector q, the integral equations 
© are transformed to a system of algebraic equations. 
Since, only a few harmonics appear in such an expansion 
of r q (k, k'), essentially those which are present in 7 q , 
the resulting algebraic system can be easily solved ana- 
lytically. We shall now describe further details for the 
two models separately. 
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FIG. 1: Graphical representation of the integral equation for 
four-point vertex. 

Model A: Frustrated antiferromagnet on a square lat- 
tice with the nearest-neighbor exchange constant J\ = 1 
and the diagonal coupling of strength J2 — \ J\ has in- 
finitely many classical ground states for < H < H c 
[H EH. The magnon spectrum in the ferromagnetic 
phase at H = H c is e q = (1 + cosg x )(l + cosq y ). The 
magnon energy has lines of minima spanned by wave- 
vectors {q*}: (tt, q) and (q, tt). Besides the well-known 
singularity related to vanishing of a scattering amplitude 
for two quantum particles in two dimensions (2d) [12| , the 
present Id type degeneracy leads to an extra infrared di- 
vergence in the kernel of the integral equation when the 



external momenta are fixed to k, k' = q* . Away from 
the critical point /i > the singularity can be cured by 
introducing an infrared cut-off for single-particle energies 
defined by £k,k' > A 4 [13 • The physics behind this reg- 
ularization procedure is as follows: at the energy scale 
of the order of /1 the interaction effects become impor- 
tant and modify the form of the bare spectrum. As it 
follows, an infrared behavior of the excitation spectrum 
is sound-like and, hence, removes in a self-consistent way 
the divergence of the kernel below H c . 

The condensate wave-vector is chosen by considering 
the scattering amplitude at zero momentum transfer and 
external momenta set to q*. The two main candidate 
states are the Neel state with Q = (tt, tt), which is sta- 
ble for a frustrated square lattice antiferromagnet with 
weaker diagonal bonds J2 < \J\, and the columnar state 
ordered on (tt,0) [(0, 7t)], which is stable for stronger di- 
agonal bonds J2 > \J\ Instead of presenting a 
rather lengthy solution of the algebraic system we ob- 
tain an approximate but physically transparent solu- 
tion by neglecting the momentum-transfer dependence 
of r p (q*,q*) = T p (q*) near p = under the integral 
in equation The result is r (q*) » l/r(q*) with 
r(q*) = 1/N J2 P <Wkq*+p + e q*-p]- I* is evident that 
ro(q*) is smallest for a q* at which the kernel r(q*) 
has the strongest divergence. This singles out the Neel 
wave- vector Q = (tt, tt), which has the softest excita- 
tions £Q+k — k 2 k 2 . Hence, for /i > magnons con- 
dense at q* = Q and a transverse antiferromagnetic 
order should be formed below the saturation field. Di- 
rect estimation of the integral in the kernel r(q*) yields 
To(Q) ~ A* 1 ^ 2 / log A* 1 1 whereas away from ordering wave- 
vector r (q*) ~ fi 1 / 2 . The square root behavior of the 
four-point vertex is intrinsic to Id systems. This is a 
consequence of the quasi Id low-energy part of the bare 
spectrum of the frustrated model. In addition, there is 
an extra logarithmic correction to the Id behavior at the 
ordering wave-vector Q related to vanishing stiffness in 
both directions leading to a Van Hove type singularity in 
the kernel. 

In the ordered state below H c , the condensate den- 
sity n = (&q) 2 is found from minimization of the 
ground state energy density e g . s . = — {ino + iro(Q)riQ: 
ri = /i/r (Q) — H 1 / 2 ] log /x I, where we have neglected 
the noncondensate contribution. The excitation spec- 
trum can be obtained following the standard Bogoliubov 
scheme and replacing the bare vertex Vq with the full 
scattering amplitude: 



q 

£ q x =n [r (q, Q) + r q _ Q (q, Q)], E q 2 =n T^(Q). 

As it is seen from Eq. (£Q) the self-energy acquires an addi- 
tional momentum dependence thanks to the dependence 
of r on the incoming momenta. The magnon spectrum 
Q is no longer degenerate and has the unique zero- 
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energy mode at q = Q. For q — > Q the magnon energy 
becomes w q « \J 4/i (q, Q) — fj]. Expanding the ver- 
tex in small |k| = |q — Q| we obtain an acoustic mode 
~ sk with the velocity s 2 ~ log/x|, which is smaller 
than the velocity s 2 ~ fi in nonfrustrated 2d and 3d an- 
tiferromagnets near the saturation. Away from the or- 
dering wave- vector, along the degeneracy lines, magnons 
acquire a dynamically generated gap A ~ u\ log/x|. 

The magnetization of the spin system is related to the 
total density of particles M(H) = h — (no + n'), which 
includes both the condensate no and noncondensate n' 
parts. The noncondensate density of bosons is given by 
n' = l/N E q [ £ q - M + Sq 1 - ^ q ]/[2w q ]. The largest con- 
tribution to n' is determined by wave- vectors away from 
Q, along the degeneracy lines q 6 {q*}- The ratio of 
the two densities is estimated as n'/n a ~ 1/llogul 1 / 2 
and is logarithmically small, though it exceeds the con- 
densate depletion found for a nondegenerate 2d Bose gas 
(n'/no ~ 1 / 1 log yu. | ) [l^ . Near H c , the magnetization 
curve exhibits a strong singularity 



M(H) ~ /x 1/2 |log/i| ~ VH C -H \\og(H c -H)\ , (5) 



which fits well to the available numerical data 

The above results obtained in the hard-core boson 
picture are valid in the vicinity of H c and should be 
compared with the standard linear spin-wave theory 
(LSWT), which applies for all fields < H < H c . In 
the LSWT approach one selects a few classical ground 
states and calculates energies of zero-point oscillations 
E n = ^X^k^k' which are different for each state. In 
this way we have found that the Neel state has lower en- 
ergy than a collinear (stripe) state in the whole range of 
fields. Another candidate state for \H C < H < H c is a 
four-sublattice state, which has identical polarization in 
three sublattices with the fourth sublattice compensat- 
ing the net transverse magnetization. Such a partially 
collinear state is a natural favorite if quantum effects 
are taken into account via effective biquadratic exchange 
interaction derived in a second-order real-space pertur- 
bation theory [l^. LSWT shows instead that the Neel 
state has again a lower energy for a frustrated square- 
lattice antiferromagnet. This discrepancy is explained 
by a nonanalytic dependence of the ground state energy 
on applied magnetic field determined by a large number 
of soft (zero) modes, whereas the real-space perturba- 
tion approach captures only analytic contributions. Fur- 
thermore, the Neel antiferromagnetic order is stable even 
at H = \H C , where it has a lower energy of zero-point 
oscillations Eq = 0.694S* than a fully collinear up-up- 
up-down (uuud) state with Eq = 0.703S*. This result 
corrects the previous spin-wave calculation an d is in 
agood agreement with the exact diagonalization studies 
|s |Uj, which show that the 1/2-magnetization plateau 
appears in the present model for 0.5 < J2/J1 ^ 0.65. 
In the vicinity of H c , the LSWT magnetization curve for 
the Neel state shows the same type of singularity (jSJ) with 



M = Sh{l + ^2(1 - h) ln[128(l - /i)/7r 4 e 2 ]/7r 2 5}, where 
h = H/H c . 

Model B: A nearest- neighbor antiferromagnet (J = 1) 
on an fee lattice is another frustrated model with de- 
generate ground states at < H < H c . The spec- 
trum of magnon excitation at the saturation field is 
e q = 2[1 + cos q x cos q y + cos q x cos q z + cos q y cos q z ] and 
has lines of zeros at (it, q, 0) and the cubic symmetry 
related lines. We then follow the same scheme as for 
the previous model and examine the kernel for the scat- 
tering amplitude with the new spectrum. The low- 
est value of renormalized four-point boson vertex corre- 
sponds to the three wave- vectors = (n, 0, 0), (0, n, 0), 
and (0, 0, 7r). The corresponding vertex is estimated as 
T = r (Q,:,Qi) ~ l/|log/i| 2 . It has a 2d logarithmic 
behavior due to a quasi-2d form of the spectrum in the 
vicinity of (n, q, 0) line with an extra logarithmic singu- 
larity related again to a vanishing stiffness for k = 0: 
£Q 1+ k ~ fc 2 + \kyk~l. Due to a presence of three equally 
singular wave-vectors one must check now whether the 
ground state of the systems is characterized by a Bose 
condensation at a single wave-vector or at all three 
wave- vectors simultaneously |l4j . For this we write the 
ground state energy in the Landau form: 



E = 



(6) 



E 



where ipi = is a complex order parameter, T = 

r (Q*, Qi) + r Q! - Qj (Q h q,-), and f = r Qi _ Q3 (q,, q,). 

A single component phase is stabilized for a sufficiently 
strong repulsion between components T < T — \T\, 
whereas in the opposite case all three components of the 
Bose condensate appear with equal weights. Direct cal- 
culations show that f ~ l/|log/i| and T ~ l/|log/i| 2 . 
Hence, in the vicinity of the saturation field dominated 
by the logarithmic behavior r » T, T and the single-k 
state is energetically favorable. The boson density and 
asymptotic behavior of the magnetization curve is given 
byn = l/2-M(i?) - //|log/i| 2 . We again find a d-l like 
form for the magnetization with a logarithmic correction. 

Corrections beyond the leading logarithms can, how- 
ever, change the energy balance. From explicit esti- 
mate of the prefactors, one finds that T > T only in 
an extremely narrow interval near the saturation field 
for /i < 7r 2 e~ 20 . Beyond this interval f < T and a 
multi-k state with a real (T < 0) superposition of all 
three modes with equal amplitudes is stabilized as the 
ground state. Such a spin structure corresponds to a 
partially collinear spin configuration described above in 
our discussion of the model A. The LSWT calculations 
confirm the above conclusion showing that for magnetic 
fields as close to the saturation field as AH/H C = 10~ 3 
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the partially collinear state is energetically more favor- 
able than a single-k state. In contrast to the behavior 
of a frustrated square-lattice antiferromagnet we have 
also found an uuud configuration for the ground state 
at H = \H C . This state has the zero-point oscilla- 
tion energy of E = 1.665, while a single-k state has 
Eq = 1.745. Thus, a 1/2-magnetization plateau should 
appear on the magnetization curve of an fee antiferro- 
magnet. 

Finally, we make a few remarks on finite temperature 
behavior. The models A and B are magnetic analogs 
of the weak-crystallization model 0]. In this model 
the phonon spectrum softens for wave-vectors lying on 
a sphere. The phase transition to a crystal state appears 
to be of the first order even if a mean-field theory would 
predict a continuous transition. Such a fluctuation driven 
first-order transition is explained by a singular Hartree 
correction generated by thermal fluctuations, which have 
a large phase space 

For a 3d quantum fee antiferromagnet a finite tem- 
perature transition to the ordered state is expected for 
all H < H c . In zero magnetic field neglecting quantum 
fluctuations the Hartree correction is estimated as Sh ~ 
E q 77[e q + r] ~ CT\ logr| 2 , where r = [T - T M f]/T M f 
is a distance from a mean- field critical point Tmf- The 
self-consistent equation for the excitation gap becomes 
A = r + CT\ log A | 2 . The point of absolute instability 
A = cannot be reached at any finite T indicating a 
first-order transition. For classical Heisenberg and XY 
model on an fee lattice the first-order transition at H = 
has been confirmed by the Monte Carlo simulations . 
The transition between paramagnetic and antiferromag- 
netic states is, therefore, of the first order at low H and 
at high T. The zero temperature phase transition at 
H c is instead of the second order. This opens two pos- 
sible scenarios for the H—T phase diagram (a) a finite 
temperature tricritical point separates low-field (high- 
T) first-order transition line from the XY transition at 
high-fields; or (b) fluctuation induced first-order transi- 
tion survives down to zero temperature and terminates at 
the quantum critical end point. For the two-dimensional 
model A, similar questions can be raised on the interplay 
between a Berezinski-Kosterlitz-Thoulees (BKT) transi- 
tion and a fluctuation induced first-order transition. At 
high temperatures (low fields) high degree of degeneracy 
of the spectra may induce the first-order transition and 
make the BKT instability inaccessible. Such a scenario 
is, for example, realized in the theory of weak crystal- 
lization of films While at low temperatures (high 
fields) a usual BKT transition may take place. These 
interesting issues require further investigations. 

Summary: We have studied a quantum phase tran- 
sition at high magnetic fields for a class of frustrated 
antiferromagnets in d = 2 and d = 3, which have de- 
generate excitation spectra with lines of minima. Mo- 



mentum dependence of the exact four-point boson vertex 
removes the degeneracy and selects the ordering wave- 
vector Q. The new spectra are sound-like near Q and 
acquire dynamically generated gaps away from it. The 
asymptotic behavior of the magnetization curve shows 
the same singularity as a nonfrustratcd model in d — 1 
dimension with an additional logarithmic correction. The 
developed scheme can be applied to other frustrated an- 
tiferromagnets near the saturation and to singlet ground 
state systems with degenerate gapped triplet excitations, 
as, e.g., SrCu 2 (B0 3 ) 2 or Cs 3 Cr 2 Br 9 0, near the 
triplet condensation field H c ±. In addition, recent exper- 
imental progress on the Bose condensation of alkali atoms 
on optical lattices |20| opens a new possibility for experi- 
mental study of condensate phenomena on frustrated lat- 
tices, the systems to which our results also apply. 
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